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INTRODUCTION 

NUSC  has  been  involved  in  the  development  of  a  Loran-C  antenna  and  preamplifier/ 
matching  network  that  is  compatible  with  a  commercial  GPS  antenna  and  navigation  receiver. 

To  be  compatible  with  the  navigation  receiver,  the  Loran  antenna  must  be  linearly  polarized  and 
provide  omni-azimuthal  reception  at  the  operational  100  kHz  carrier  frequency  with  a  20  kHz 
bandwidth.  In  addition,  the  Loran  antenna  is  electrically  small  since  it  is  designed  to  fit  within  a 
cylindrical  volume  of  maximum  height  between  4  and  8  in.  and  maximum  diameter  between  3 
and  8  in.  The  restrictions  on  the  maximum  dimensions  of  the  antenna  were  set  during  the 
analysis  stage  of  the  design.  The  Loran-C/GPS  package  is  designed  to  be  placed  on  a  floating 
platform.  Fig.  1  gives  a  block  diagram  of  the  Loran-C/GPS  integrated  antenna  system. 

This  report  is  primarily  concerned  with  the  theoretical  determination  of  an  optimum  antenna 
for  Loran-C  reception  that  fits  within  the  allowable  dimensions.  This  study  includes  an  analysis 
of  the  amplifier  noise  along  with  a  determination  of  the  input  reactance  and  the  effective  height  of 
the  antenna.  The  analysis  is  restricted  to  monopole  antennas  since  a  constant  amplitude  and 
phase  response  is  required  for  all  azimuthal  angles.  Wherever  possible,  the  theoretical  models  are 
validated  by  experimental  data. 


CIRCUIT  ANALYSIS 

Consider  an  electrically  short  monopole  of  height  h  above  a  ground  plane  as  shown  in  Fig.  2a. 
The  Thevenin  equivalent  circuit  of  the  monopole  receiving  a  signal  is  given  in  Fig.  2b,  where 
VQC  =  open  circuit  voltage  received  by  the  antenna, 

Ra  =  radiation  resistance  of  the  antenna, 

R}  =  ohmic  resistance  of  the  antenna, 

Ca  =  antenna  capacitance. 

Rj  includes  the  ohmic  losses  associated  with  the  finite  conductivities  of  the  antenna  and  the  sea 
water  ground  plane.  Since  the  loss  tangent  in  sea  water  is  very  large  over  the  frequency  range  of 
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concern,  the  ground  plane  may  be  approximated  by  a  perfect  conductor.  The  input  impedance  of  the 
monopole  antenna  is  R  +  l/jcoCa,  where  to  is  the  angular  frequency  while  R  =  Ra  +  Rj  and  R  « 
l/j(oCa.  The  antenna  capacitance  is  obtained  by  solving  an  electrostatics  problem  [1].  In  the 
limiting  case  of  a  vanishingly  small  frequency,  a  dipole  or  monopole  antenna  is  just  a  capacitor. 
For  a  uniform  incident  electric  field  Emc  directed  parallel  to  and  in  the  vicinity  of  the  antenna,  the 
open  circuit  voltage  is  given  by  [2] 

(1) 

where  he^  is  the  effective  height  of  the  antenna.  The  computation  of  the  antenna  capacitance  and 
effective  height  will  be  discussed  later. 

Next  consider  the  monopole  connected  by  a  short  wire  to  the  input  of  an  FET  amplifier  as  shown 
in  Fig.  3a  with  equivalent  circuit  given  in  Fig.  3b,  where 

Cw  =  interconnection  capacitance  between  the  antenna  and  the  amplifier, 

R0  =  input  resistance  of  the  amplifier, 

CQ  =  input  capacitance  of  the  amplifier, 

Vn  =  open  circuit  noise  voltage  produced  by  .the  amplifier 
To  compare  the  total  noise  produced  by  the  amplifier  with  the  atmospheric  noise  received  by  the 
antenna,  all  voltages  must  be  referred  to  the  same  location.  For  convenience,  we  will  refer  all 

noise  voltages  to  the  antenna  terminals.  The  thermal  noise  voltage  associated  with  the  loss 
resistance  Rj  of  the  antenna  has  been  computed  to  be  of  negligible  value  and  will  not  be  considered 

further  in  this  analysis. 

To  determine  the  amplifier  noise,  consider  Figs.  4(a)  and  4(b).  The  voltage  En  and  current  In 

are  equivalent  noise  sources  that  adequately  represent  the  noise  produced  by  the  amplifier  and 

whose  values  can  easily  be  determined  experimentally  [3].  The  total  noise  power  produced  by  the 
amplifier  is  given  by  the  sum  of  the  contributions  from  En  and  In.  Vjj  is  the  voltage  produced  by  En 

at  the  input  terminals  of  the  amplifier  while  Vi2  is  the  voltage  produced  by  In  at  the  same  location. 
Ya>  Yw>  and  Y0  are  the  admittances  of  the  antenna,  the  interconnection,  and  the  amplifier, 
respectively.  An  analysis  of  Figs.  4(a)  and  4(b)  yields 
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Vii  =  En 


Ya+Yw 
Ya+ Yw  + Y0 


(2a) 


Vi2  =  In 


Ya+Yw  +  YJ 


(2b) 


Vii  and  Vi2  can  be  referred  to  equivalent  sources  at  the  antenna  terminals  by  considering  Fig.  5. 
In  Fig.  5,  a  voltage  source  Vs  at  the  antenna  terminals  produces  a  voltage  Vig  at  the  input  of  the 
amplifier.  An  analysis  of  Fig.  5  yields 


V-  Y 

Tr  is _ a 

Y  ~  Y  +  Y  +  Y 

s  a  w T  o 


where  Ky  is  the  voltage  transfer  function.  The  application  of  (3)  to  (2a)  and  (2b)  gives  the 
equivalent  amplifier  noise  voltage  sources  Vnl  and  Vn2  due  to  En  and  In,  respectively,  at  the 
antenna  terminals.  Therefore,  from  (2)  and  (3), 


(3) 


(4a) 


v  -Zii-T  J_ 

02  Kv  n  Ya 

If  Vnl  and  Vn2  are  uncorrelated,  the  open-circuit  noise  voltage  Vn  produced  by  the  amplifier  is 
given  by 

172 


(4b) 


7  -  h 

\l 


v„  =  iivnl|2+|vn2|2 


,2 

,  Y  w 

2 

It  |2 

1 

2 

eJ 

1  + 

+  In 

ni 

Ya 

1  "I 

Ya 

,1/2 


(5) 


Note  that  Vn  is  independent  of  the  amplifier  impedance. 

For  convenience,  we  may  relate  the  amplifier  noise  voltage  to  an  equivalent  incident  electric 
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field  through  the  effective  height  of  the  receiving  antenna.  From  (1)  and  (5), 


amp  Vn 


(6) 


gamp  incident  electric  field  parallel  to  the  antenna  that  produces  the  same  noise  power  as  the 
amplifier.  The  substitution  of  Ya  =  jcoCa  and  Yw  =  jcoCw  into  (6)  gives 


E^p 


(7) 


The  goal  of  this  effort  is  to  design  an  antenna  that  produces  a  minimum  value  of  Eamp.  Therefore, 
it  is  desirable  to  maximize  both  the  antenna  capacitance  and  the  effective  height.  This 
requirement  results  in  a  tradeoff  that  will  be  discussed  later.  It  is  preferred  that  the 
antenna/amplifier  system  be  atmospheric  noise  limited;  i.e.,  Eamp  <  Eatm,  where  Eatm  is  the 
atmospheric  noise  field  as  determined  from  the  CCIR  manual  [4]. 


ANTENNA  CAPACITANCE 

The  electrostatic  charge  distribution  and  capacitance  of  an  electrically  small  antenna  play  an 
important  role  in  determining  its  circuit  and  radiation  characteristics.  No  analytical  expression 
exists  for  the  charge  distribution  along  an  arbitrary  conducting  body  of  finite  dimensions  in  free 
space  or  over  a  ground  plane.  This  section  presents  a  numerical  algorithm  for  determining  the 
electrostatic  charge  density  distribution  and  capacitance  of  a  conducting  body  of  revolution  above 
a  perfectly  conducting  ground  plane.  The  algorithm  is  applied  to  the  analysis  of  several  monopole 
antenna  geometries  including  cylinders  and  truncated  cones.  Each  structure  will  be  studied  both 
with  and  without  a  top  load.  The  algorithm  is  verified  by  comparing  the  computed  capacitances  of 
various  cylinders  with  measured  data. 

INTEGRAL  EQUATION  FORMULATION 

Consider  a  body  of  revolution  which  is  formed  by  rotating  a  planar  curve  about  the  z-axis  as 
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described  in  Pig.  6.  In  Fig.  6,  (p,<|>,z)  are  cylindrical  coordinates,  and  (t,c(>)  form  an  orthogonal 
curvilinear  coordinate  system  on  the  surface  S  of  the  body  of  revolution.  The  coordinate  t  is  the  arc- 
length  along  the  generating  curve  and  <{>  is  the  azimuthal  angle  measured  from  the  x-z  plane.  The 
ends  of  the  generating  curve  may  touch  the  z-axis. 

If  the  body  of  revolution  is  a  conductor  charged  to  a  potential  V  with  respect  to  a  ground  plane 
located  on  the  x-y  plane,  the  following  integral  equation  for  the  unknown  surface  charge  density 
cts  results  [1]: 

r  r 


i 

4rte0 


Jj  ffrfr') 


dS'  =  V 


re  S 


(8) 


where  r  and  r'  are  position  vectors  of  the  observation  and  source  points,  respectively,  and  the 
integration  is  over  the  surface  of  the  body.  In  (8),  R+  and  R-  denote  the  distances  between  the 
observation  and  source  points  on  the  body  and  its  image,  respectively.  Since  the  boundary 
conditions  and  excitation  are  axisymmetric,  all  field  quantities  are  (^-independent.  Therefore,  <|> 
may  be  set  equal  to  zero  and  suppressed.  In  the  coordinate  system  of  the  surface  of  revolution,  (8) 
may  be  expressed  as 


p’(t')  d<t>'  dt'  =  V 


reS 


(9) 


where 


dt'  =  [  (dp')2  +  (dz’)2  ]m 


(10) 


R±=  Ir+r’l  =  [  (p-p’)2  +  (z  +  z')2  +  4pp'  sin2|-] 


1/2 


(ID 


T  is  the  total  am  length  along  the  generating  curve.  The  unprimed  coordinates  refer  to  the 
observation  point  while  the  primed  coordinates  refer  to  the  source  points.  Note  that  p  and  z  can  be 
expressed  as  functions  of  t  where  (p,z)  and  (p',z')  are  the  cylindrical  coordinates  of  the  points  t  and 
t',  respectively. 
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Since  the  surface  charge  density  is  axisymmetric,  we  may  write 

q(t')  =  27ip(t')crs(t') 

where  q(t)  is  the  linear  charge  density.  The  substitution  of  (12)  into  (9)  gives 


(12) 


_L 


4jt2e0 


fT 

"  f* 

q(t') 

Jo 

]0  (£-x)df 

dt’  =  V 


(13) 


Note  that  the  even  symmetry  of  the  integrand  with  respect  to  4»'  =  0  has  been  accounted  for  in  (13). 
1 

The  substitution  a  =  ^(n-<t>')  in  the  <t»'  integration  of  (13)  yields 


f 


R+  R 


dp1  =  -2-  K(u2)  --2-  K(w2) 

R1  R7 


where 


Rf  =  [  (p  +  p’)a  +  (z  +  z?  ] 
2^PP7 


u  = 


R1 


w  - 


(14) 

(15) 

(16a) 

(16b) 


kl 


en/2 

K(m)  =  I  d« 

Jo  Vl- 


(17) 


.  2 

m  sm  a 


K(m)  is  the  complete  elliptic  integral  of  the  first  kind  [5].  Hence,  from  (14),  (13)  can  be  rewritten  as 

/■T 


2 

2n  t, 


q(t') 


K(u2)  K(w2) 


R, 


R, 


dt1  =  V 


(18) 


Expression  (18)  is  an  integral  equation  for  the  unknown  linear  charge  density  q(t)  along  the 
conducting  body  of  revolution.  The  advantage  of  (18)  over  (13)  is  that  there  exist  polynomial 
approximations  for  the  elliptic  integrals  which  can  be  integrated  with  ease  [51.  Once  q(t)  has  been 
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determined,  the  capacitance  C  of  the  body  with  respect  to  ground  can  be  evaluated  as 

•  T 


_  Q  =!f 
v'vl 


q(t')dt' 


(19) 


NUMERICAL  FORMULATION 

Consider  the  approximation  of  the  generating  curve  of  the  body  of  revolution  as  a  sequence  of  N 

linear  segments  as  shown  in  Fig.  7.  The  endpoints  of  the  segments  lie  on  the  generating  curve  and 
are  denoted  by  tg,  tj,  t2,  tjj.  Each  endpoint  is  described  in  cylindrical  coordinates  as  t-  =  (pj.zj) 

i  =  0,  1,  2, N.  The  segment  midpoints  are  given  by 


**  =  2 
while  the  segment  lengths  are 


Vi +  \  /  Pi-1 +  Pi  Vi  +  zi 


=  ( 


)  =  (Pj.Zj)  i  =  1,  2, N  (20) 


Aj  =  [(pi  ~Pi-l)2  +  (zi  "Vl)2] 


V2 


i  =  1,  2, N  (21) 


On  segment  j  (  tj_j  <  t  <  tj  )  the  following  parameterization  is  valid: 

z  =  Zj_x  +  (t  -  tj.j)  cos  Vj 


P  =  pj_1  +  (t-tj_1)sinvj 


where 


Vj  =  tan 


-1 


/-  - 


PHii-ll 


(22a) 


(22b) 


(23) 


Here,  Vj  is  the  angle  between  the  approximate  generating  curve  at  t  =  tj  and  the  z-axis. 

The  method  of  moments  [6]  may  be  applied  to  obtain  a  numerical  solution  to  the  integral 
equation  (18).  In  this  implementation,  the  linear  charge  density  is  approximated  as  a  linear  sum 
of  pulses: 


N 


q(t)  =  ]T  Qj  P/t) 
j=l 


(24) 


where 
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1,  Vj  <  t  ^  tj 
0,  elsewhere 


(25) 


Qj  is  the  linear  charge  density  at  the  center  of  the  jth  segment  (t  =  tj).  Fig.  8  gives  a  depiction  of  the 


pulse  approximation  of  q(t).  The  substitution  of  (24)  and  (25)  into  (18)  followed  by  an  enforcement 
of  the  resulting  equation  at  the  center  of  each  segment  (t  =  tj,  i  =  1,  2, ...,  N)  yields  a  linear  system  of 
N  equations  in  the  N  unknown  coefficients  Qj: 


N 

^ZijQj  =  Vi  i  =  1)  2, N  (26) 

j-1 


where 


Vi  =  2*V 


K(wfj) 


R; 


dx' 


Rj  =  [  (pj  +  Pj_i  +  <  sin  Vj)2  +  (zj  +  Zj_1  -  x'  cos  Vj)2  ] 


1/2 


(27) 


(28) 


(29) 


.£ 


172 


(30a) 


_2[Pi(Pj-l  +  T’sinvj)l 

Wij  "  - - - - 

Rij 


(30b) 


The  integral  in  (28)  was  obtained  following  the  change  of  variables  x'  =  t'  -  tj_j.  In  (28),  the  first 

term  of  the  integrand  possesses  a  logarithmic  singularity  when  the  observation  point  is  within  the 
source  region  (i  =  j).  From  Appendix  A, 


K(u2) 


-^-[ln(4)  +  ln(2p)-lnlt-t'l] 


(31) 
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Since  only  the  last  term  in  (31)  is  singular,  we  may  subtract  it  from  the  integrand  and  add  its 
analytically  evaluated  integral.  After  an  appropiate  change  of  variables,  we  have  from  (28)  and 


(31), 


Ai 


Zji  = 


K(ufj)  K(wJ)  1 


R; 


Ri 


+  —  ln|  — -x’ 

2Pi  ‘ 


-irfH*- 


dx’  (32) 


The  second  integral  in  (32)  may  be  evaluated  analytically  as 

,Ai 


2p 

The  substitution  of  (33)  into  (32)  gives 

Ai 


Zii  = 


dx'=-l. 

2P1 


lnl—  -1 


K<"ii>  K(wii)  ,  1  IlJAi  T’ 

R+.  R-  2pi  I  2 


dx'  +  — 
2Pi 


1-  ln|— 
2 


The  representation  of  as  given  by  (34)  can  be  integrated  numerically. 

The  linear  system  of  equations  (26)  may  be  written  in  matrix  form  as  follows: 

[Z][Q]  =  [V] 


(33) 


(34) 


(35) 


where  [Z]isanNxN  coefficient  matrix  with  elements  given  by  (28)  and  (34)  while  [  V  ]  is  the 
excitation  vector  with  elements  given  by  (27).  The  unknown  charge  density  vector  [  Q  ]  =  [  Qj,  Q2, 

T 

Qn  1  is  determined  by  matrix  inversion.  An  explanation  and  listing  of  a  FORTRAN 


language  computer  program  based  on  the  method  of  moments  solution  of  (18)  as  described  in  this 
section  is  given  in  Appendix  B. 


EXPERIMENTAL  SETUP 

The  numerical  algorithm  was  tested  by  comparing  its  results  with  the  measured  capacitances 
of  several  hollow  cylinders  at  varying  distances  above  the  ground  plane.  Various  tubes  were 
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constructed  with  heights  varying  from  1  in.  to  7  in.  and  diameters  ranging  from  3.5  in.  to  8.2  in. 
The  capacitances  were  measured  at  100  kHz  using  a  low  frequency  network  analyzer.  The  ground 
plane  used  in  this  investigation  was  a  3  ft  diameter  sheet  of  1/4  in.  thick  steel.  These  dimensions 
were  found  to  adequately  approximate  an  infinite  ground  plane.  The  tubes  were  constructed  by 
attaching  thin  copper  tape  around  PVC  pipes  of  various  diameters.  A  small  piece  of  wire  was 
soldered  to  the  outside  of  each  tube  for  connection  of  the  network  analyzer  probe.  The  distance 
above  the  ground  plane  was  varied  by  placing  the  tube  on  top  of  several  small  wooden  blocks  of 
different  thicknesses.  The  entire  measurement  setup  was  placed  in  a  shielded  room. 

The  tubes  were  also  measured  with  a  top  load.  The  top  load  consisted  of  an  8  1/4  in.  diameter 
flat  disk  of  copper  tape  attached  to  a  cardboard  backing  that  was  soldered  to  the  tube  under  study. 
The  measurement  repeatability  was  found  to  be  approximately  ±.75  pF  for  the  unloaded  tubes  and 
±1.5  pF  for  the  top-loaded  tubes. 

RESULTS 

Consider  a  highly  conducting  tube  of  length  h,  radius  a,  and  separated  by  a  distance  s  above  a 

perfectly  conducting  ground  plane  as  illustrated  in  Fig.  9(a).  Fig.  10  shows  the  normalized  linear 
charge  density  q(t)/e0V  of  a  tube  of  length  h  =  6  in.  and  diameter  2a  =  3.5  in.  plotted  as  a  function  of 

axial  distance  t  along  the  tube  for  various  separations  above  the  ground  plane.  The  curves  show  a 
charge  singularity  at  each  end  of  the  tube  as  predicted  from  theory  [7],  For  large  separations  above 
the  ground  plane,  the  charge  distribution  is  symmetric  about  the  midpoint  of  the  tube.  As  the  tube 
separation  above  the  ground  plane  becomes  small,  the  charge  symmetry  disappears  with  more 
charge  accumulating  on  the  ground  plane  side  of  the  tube. 

Fig.  11  shows  a  comparison  of  the  computed  and  measured  capacitances  of  several  6  in.  long 
conducting  tubes  of  different  radii  plotted  as  a  function  of  separation  above  the  ground  plane.  The 
curves  show  that  the  capacitance  decreases  with  separation  from  the  ground  plane  and  increases 
for  larger  radii.  The  data  indicate  that  the  measured  and  computed  capacitances  agree  to  within 
1  pF  for  most  separations,  an  excellent  correspondence  considering  the  margin  for  error  in  the 
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measurements.  A  similar  comparison  of  the  computed  and  measured  capacitances  of  several  3  in. 
long  tubes  given  in  Fig.  12  also  shows  an  agreement  to  within  1  pF.  Comparisons  for  tubes  of  other 
heights  produced  similar  results,  verifying  the  accuracy  of  the  body  of  revolution  algorithm. 

Appendix  C  gives  an  approximate  formula  (C  -  3)  for  half  of  the  capacitance  of  an  electrically 
small  tubular  monopole  above  a  perfectly  conducting  ground  plane.  This  formula  is  a  modified 
expression  for  the  capacitance  of  a  coplanar  stripline.  The  advantage  of  this  formula  over  the 
moment-method  solution  is  that  its  evaluation  involves  only  the  computation  of  elliptic  integrals 
of  the  first  kind  for  which  simple  polynomial  approximations  exist  [5].  A  comparison  of  the 
approximate  formula  with  the  integral  equation  solution  and  measured  data  for  the  capacitances 
of  several  5  in.  long  conducting  tubes  is  given  in  Fig.  13.  The  curves  show  a  close  agreement  with 
both  the  measured  and  integral  equation  data,  indicating  the  adequacy  of  the  modified  coplanar 
stripline  formula  for  computing  the  capacitance  of  an  electrically  small  monopole  or  dipole 
antenna.  For  additional  results,  the  reader  is  referred  to  [8]. 

Consider  the  top-loaded  tube  shown  in  Fig.  9(b).  It  consists  of  a  highly  conducting  thin  disk  of 
radius  b  that  is  attached  to  the  tube  in  Fig.  9(a).  Since  the  geometry  of  the  top-loaded  tube  is  not 
tractable  using  the  body  of  revolution  algorithm  described  in  the  previous  section,  the  top  load  and 
the  tube  had  to  be  analyzed  as  two  coupled  bodies  of  revolution.  The  specifics  of  a  revised  algorithm 
to  analyze  this  structure  will  not  be  given  in  this  report  as  it  is  a  direct  extension  of  the  previously 
described  body  of  revolution  algorithm. 

Fig.  14  shows  a  comparison  of  the  computed  and  measured  capacitances  of  several  top-loaded 
conducting  tubes  with  h  =  6  in.  and  2b  =  8.25  in.  of  different  radii  plotted  as  a  function  of  separation 
above  the  ground  plane.  The  data  show  the  same  behavior  as  the  unloaded  tubes  in  Fig.  11  except 
that  the  capacitances  are  larger.  The  data  also  show  more  deviation  between  the  measured  and 
computed  results.  This  is  expected  since  the  measurement  repeatability  is  worse  in  this  case.  A 
similar  comparison  of  the  computed  and  measured  capacitances  of  several  top-loaded  tubes  with  h 
=  3  in.  and  2b  =  8.25  in.  is  given  in  Fig.  15.  Comparisons  for  top-loaded  tubes  of  other  heights  with 
2b  =  8.25  in.  produced  similar  results,  further  verifying  the  accuracy  of  the  revised  body  of 
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revolution  algorithm. 

The  body  of  revolution  algorithm  was  also  applied  to  compute  the  static  charge  distribution  and 
capacitance  of  a  conical  antenna  over  a  conducting  ground  plane.  Our  computed  data  compare 
well  with  Wilton's  computed  data  [9]  for  various  conical  antennas  both  with  and  without  a  top  load. 

When  the  body  of  revolution  described  in  Fig.  6  is  used  as  an  antenna,  it  must  be  connected  to  a 
receiver  located  below  the  ground  plane.  Therefore,  the  portion  of  the  connecting  wire  above  the 
ground  plane  and  the  body  of  revolution  form  an  antenna  which  acts  as  a  monopole.  Since  the 
connecting  wire  is  at  the  same  potential  as  the  body  of  revolution,  the  antenna  capacitance  is  the 
sum  of  the  wire  and  body  capacitances.  If  the  connecting  wire  is  very  thin  such  that  it  carries  a 
negligible  charge  as  compared  with  the  total  charge  on  the  body,  its  capacitance  may  be  ignored.  In 
this  report,  the  capacitance  of  the  connecting  wire  will  be  neglected. 


EFFECTIVE  HEIGHT 


The  effective  height  heg-  of  a  linearly  polarized  antenna  receiving  a  plane  wave  from  a  given 

direction  is  the  ratio  of  the  magnitude  of  the  open  circuit  voltage  developed  at  the  terminals  of  the 
antenna  to  the  magnitude  of  the  electric  field  strength  in  the  direction  of  the  antenna  polarization 
[10].  This  definition  has  been  applied  in  (1).  Alternatively,  the  effective  height  is  the  length  of  a 
thin  straight  conductor  oriented  perpendicular  to  the  given  direction  and  parallel  to  the  antenna 
polarization,  having  a  uniform  current  equal  to  that  at  the  antenna  terminals  and  producing  the 
same  far-field  strength  as  the  antenna  in  that  direction  [10].  The  application  of  this  latter 
definition  to  the  body  of  revolution  in  Fig.  6  yields 


I(t)  cosv(t)  dt 


(36) 


where  T  is  the  arc  length  along  the  generating  curve  of  the  antenna,  v(t)  is  the  angle  between  the 
tangent  to  the  generating  curve  and  the  z-axis,  I(t)  is  the  current  along  the  antenna,  and  1(0)  is  the 
current  at  the  antenna  terminals.  It  is  assumed  that  the  body  is  separated  by  an  infinitesimal 
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distance  from  the  ground  plane. 

An  expression  for  the  effective  height  of  a  body  of  revolution  may  be  obtained  in  terms  of  the 
charge  distribution.  The  application  of  the  method  of  moments  formulation  described  in  the 
previous  section  to  the  effective  height  formula  (36)  gives 


From  the  continuity  equation, 

,  I'(t)  =  -jci)q(t) 

The  integration  of  (38)  over  the  entire  antenna  yields 


I(T)-I(0) 


r 

=  -j®  | 

«f\ 


q(t)dt 


(37) 


(38) 


Since  I(T)  =  0, 


T  •  N 

1(0)  =  jco  I  q(t)dt  =  jco  V  QjAj 

i=l 


(39) 


where  Qj  A-  is  the  total  charge  on  segment  i.  The  application  of  (38)  to  the  integral  in  (37)  gives 

■*i 


f 

Jk-i 


I(t)  dt  =  tiKtj)  -  tj_i  Ktj.i)  +  Qi  Ai  (ti-i + tj  ( 40 ) 


Upon  repeated  application  of  (38), 


N 


I(tj)  =  j®  ^  QjAj 
j=i+i 


(41) 


The  substitution  of  (41)  into  (40)  yields 


N 

I(t)  dt  =  j(0  Aj^  YijQjAj 

j  =i 


where 


(42) 
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r..J  m'  i=j 

>j  l  1,  otherwise 


The  substitution  of  (39)  and  (42)  into  (37)  gives 


heff  = 


N  N 

XAicosviSYiiQiAi 

i=i  j=i 


N 

i*  1 


(43) 


N  i 

XQiAi  X  YijAJC0SVj 

i  =  l  j  =  l 


N 


£  QiAi(zi~zo) 
i  =  l 


N 

is  1 


N 


(44) 


Equation  (44)  indicates  that  the  effective  height  of  an  electrically  small  body  of  revolution  is 
actually  the  axial  height  of  the  center  of  charge.  The  above  expression  also  shows  the  importance  of 
capacitive  loading  at  the  top  of  the  antenna  since  it  raises  the  center  of  charge  and  thus  increases 
the  effective  height  of  the  antenna.  A  good  discussion  of  the  effective  height  of  electrically  small 
linear  antennas,  is  given  by  Schelkunoff  and  Friis  [1]. 

The  above  formulation  applies  to  a  body  of  revolution  that  is  separated  by  an  infinitesimal 
distance  from  the  ground  plane.  If  the  body  is  separated  by  a  finite  distance  s  ( =  Zq)  above  the 

ground  plane  and  connected  to  a  receiver  (located  below  the  ground  plane)  by  a  thin  wire,  the 
effective  height  formula  (36)  must  be  modified.  If  the  connecting  wire  is  very  thin  such  that  it 
carries  a  negligible  charge  as  compared  with  the  total  charge  on  the  body  of  revolution,  the  center  of 
charge  of  the  antenna  (connecting  wire  +  body  of  revolution)  is  approximately  the  same  as  that  of 
the  body.  Hence,  the  effective  height  formula  (36)  becomes 


eff 


-mi: 


I(t)  cosv(t)  dt 


(45) 


while  the  numerical  formula  (44)  is 
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N 

XQiAi(zi“s) 

i  =  1 

heff  =  s  + - - -  <46> 

S^iAi 

i  =1 


EXAMPLES 

Fig.  16  shows  the  effective  height  plotted  as  a  function  of  diameter  for  a  conducting  tube  of 
height  h  =  6  in.  and  separated  by  0.1  in.  above  a  ground  plane.  In  this  and  all  remaining  plots,  it  is 
assumed  that  the  antenna  is  connected  to  a  receiver  (located  below  the  ground  plane)  by  an 
infinitesimally  thin  wire  so  that  (45)  and  (46)  are  valid.  As  the  tube  diameter  approaches  zero,  the 
charge  distribution  becomes  nearly  symmetrical  about  its  center,  resulting  in  an  effective  height 
of  almost  3.1  in.  As  the  tube  diameter  increases,  the  added  capacitive  coupling  from  the  lower  end 
of  the  tube  to  ground  causes  the  effective  height  to  decrease. 

Fig.  17  illustrates  a  conducting  tube  of  height  h  and  diameter  2a  that  is  top  loaded  with  a  thin 
annular  disk  of  inner  diameter  2a  and  outer  diameter  2b.  The  tube  is  separated  by  a  short  distance 
s  above  a  ground  plane.  Fig.  18  shows  the  effective  height  of  a  top-loaded  tube  with  h  =  6  in.,  2a  =  0.5 
in.,  and  s  =  0.1  in.  that  is  plotted  as  a  function  of  the  top  load  outer  diameter.  The  curve  illustrates 
the  improvement  in  effective  height  resulting  from  an  increase  in  the  size  of  the  top  load. 

Consider  the  truncated  cone  described  in  Fig.  19.  It  consists  of  a  section  of  hollow  cone  with 
axial  length  h,  lower  diameter  dj,  and  upper  diameter  df.  The  cone  is  separated  by  a  small  gap  s 

above  a  perfectly  conducting  ground  plane  and  connected  to  a  receiver  below  the  ground  plane  by 
an  infinitesimally  thin  wire.  Fig.  20  shows  the  effective  height  of  a  truncated  cone  with  h  =  6  in.,  df 
=  5  in.,  and  s  =  0.1  in.  that  is  plotted  as  a  function  of  the  lower  diameter  dj.  A  narrower  initial 
diameter  reduces  the  capacitive  coupling  from  the  lower  end  of  the  cone  to  ground,  thereby  raising 
the  center  of  charge  of  the  antenna.  Consequently,  the  curve  shows  that  the  effective  height 
increases  with  a  narrower  initial  diameter.  However,  as  will  be  shown  in  the  next  section,  this 
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increase  in  effective  height  is  attained  at  the  expense  of  a  decrease  in  capacitance. 

DETERMINATION  OF  AN  OPTIMUM  ANTENNA 

During  this  investigation,  certain  space  restrictions  for  the  Loran-C  antenna  were  set.  It  was 
decided  that  the  antenna  must  fit  within  a  cylinder  of  height  4.25  in.  and  diameter  5  in.,  where  the 
bottom  of  the  cylinder  is  coincident  with  the  sea  water  line.  The  Loran  antenna  may  be  top  loaded 
but  must  allow  for  a  feed  cable  to  pass  through  its  middle  for  connection  to  a  GPS  antenna.  Thus,  it 
was  decided  to  reserve  a  cylinder  of  diameter  0.5  in.  in  the  middle  of  the  Loran  antenna  for  the 
passage  of  the  GPS  antenna  feed.  The  space  restrictions  for  the  Loran-C  antenna  are  illustrated  in 
Fig.  21. 

Measurements  of  the  FET  amplifier  were  made  at  100  kHz  to  determine  the  equivalent  noise 
sources  En  and  In  as  defined  in  Fig.  4.  The  results  of  these  measurements  indicated  that  for 

antenna  capacitances  below  100  pF  the  amplifier  is  current  noise  limited;  i.e.,  the  noise  produced 
due  to  In  is  much  greater  than  that  due  to  En.  Therefore,  in  any  calculations  of  the  amplifier  noise 
field  as  determined  by  (7),  the  contribution  due  to  En  will  be  neglected.  The  measurements,  based 
on  a  100  Hz  bandwidth  signal  with  a  100  kHz  carrier  frequency,  indicate  that  In  is  approximately  2 

pA. 

In  this  section,  an  optimum  monopole  antenna  configuration  will  be  determined.  The 
optimum  antenna  will  be  one  that  fits  within  the  space  restrictions  as  described  in  Fig.  21  and  will 
result  in  a  minimum  amplifier  noise  field  as  computed  by  (7).  The  antenna  geometries  to  be 
considered  are  cylinders  and  truncated  cones.  Each  structure  will  be  studied  both  with  and  without 
a  top  load. 

CYLINDER 

Consider  a  tubular  monopole  as  described  in  Fig  9(a)  with  h  =  4.15  in.  and  s  =  0.1  in.  The 
antenna  capacitance,  effective  height,  and  amplifier  noise  field  are  plotted  as  functions  of  the  tube 
diameter  in  Figs.  22a,  22b,  and  22c,  respectively.  The  curves  show  that  the  antenna  capacitance  is 
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a  monotonically  increasing  function  of  the  tube  diameter  while  the  effective  height  is  a  decreasing 
function.  The  amplifier  noise  field  Eamp  as  computed  from  (7)  is  inversely  proportional  to  the 
product  of  the  effective  height  and  the  antenna  capacitance  for  a  current  noise  limited  system.  Fig. 
22c  shows  that  E8171?  is  a  monotonically  decreasing  function  of  the  tube  diameter.  Therefore,  the 
change  in  capacitance  due  to  a  variation  in  the  tube  diameter  overrides  the  opposite  change  in  the 
effective  height.  For  example,  a  tubular  monopole  with  2a  =  0.5  in.  yields  Ca  =  3.52  pf,  heff  =  1.80 
in.,  and  Eamp  =  19.8  pV/m  while  a  tube  with  2a  =  5  in.  gives  Ca  =  22.4  pf,  heg-=  1.32  in.,  and  Eamp  = 

4.24  pV/m.  Upon  increasing  the  tube  diameter  from  0.5  in.  to  5  in.,  the  capacitance  increases  by 
more  than  a  factor  of  6  while  the  effective  height  decreases  by  less  than  1.5.  Consequently,  the 
amplifier  noise  field  decreases  by  more  than  a  factor  of  4.5,  corresponding  to  an  improvement  of 
13.4  dB.  Thus,  the  data  in  Fig.  22c  indicate  that  the  optimum  tubular  monopole  is  one  with  the 
maximum  allowable  diameter  (2a  =  5  in.). 

If  the  tube  is  top  loaded  with  a  conducting  annular  disk  of  inner  diameter  0.5  in.  and  outer 

diameter  5  in.  as  described  in  Fig.  17,  some  improvement  in  performance  can  be  obtained.  For 
example,  a  top-loaded  tube  with  h  =  4.15  in.,  s  =  0.1  in.,  and  2a  =  0.5  in.  results  in  Ca  =  7.36  pf,  heff  = 
3.23  in.,  and  E8771?  =  5.27  pV/m  while  a  tube  with  2a  =  5  in.  gives  Ca  =  22.7jpf,  heff  =  1.39  in.,  and 

gamp  _  3  93  j^v/m.  Thus,  top  loading  provides  an  11.5  dB  decrease  in  the  amplifier  noise  field  for 
the  0.5  in.  diameter  tube  while  only  resulting  in  a  0.59  dB  improvement  for  the  5.0  in  diameter 
tube.  However,  the  larger  tube  still  gives  the  best  performance. 

TRUNCATED  CONE 

Consider  the  top-loaded  truncated  cone  described  in  Fig.  23.  The  antenna  capacitance, 
effective  height,  and  amplifier  noise  field  are  plotted  as  functions  of  the  lower  diameter  dj  in  Figs. 

24a,  24b,  and  24c,  respectively,  for  the  truncated  cone  both  with  and  without  the  top  load.  The  curves 
show  that  the  antenna  capacitance  increases  with  di  while  the  effective  height  exhibits  the  opposite 
behavior.  In  addition,  Fig.  24c  indicates  that  Eamp  is  a  monotonically  decreasing  function  of  dj. 
Consequently,  the  increase  in  capacitance  associated  with  an  increase  in  the  lower  diameter 
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overrides  the  corresponding  decrease  in  the  effective  height.  Therefore,  tapering  the  antenna  at  the 
base  does  not  provide  an  improvement  in  the  overall  antenna  performance.  For  example,  an 
unloaded  cone  with  dj  =  0.5  in.  results  in  Eamp  =  5.02  |iV/m  while  one  with  dj  =  5.0  in.  gives  Eamp  = 

4.24  jiV/m.  This  increase  in  the  lower  diameter  only  results  in  a  1.47  dB  improvement  in  the 
amplifier  noise  field.  Note  that  a  truncated  cone  with  dj  =  df  is  just  a  cylinder.  The  amplifier  noise 

curves  in  Fig.  24c  show  that  the  cylindrical  monopole  with  2a  =  5  in.  provides  a  lower  amplifier 
noise  field  than  a  truncated  cone. 

Fig.  24c  shows  that  a  top-loaded  cone  provides  a  lower  amplifier  noise  field  than  an  unloaded 
cone.  However,  this  improvement  is  quite  small.  For  example,  a  top-loaded  truncated  cone  with  dj 

=  0.5  in.  results  in  only  a  0.67  dB  improvement  over  an  unloaded  one. 

FINAL  DESIGN 

The  analysis  of  the  previous  section  has  shown  that  the  optimum  Loran-C  antenna  that  fits 
within  the  space  restrictions  described  in  Fig.  21  is  a  tubular  monopole  with  the  maximum 
allowable  diameter  (2a  =  5  in.).  From  Figs.  22a-22c,  this  antenna  results  in  Ca  =  22.4  pF,  he^=  1.32 

in.,  and  Eamp  =  4.24  pV/m.  Although  the  analysis  showed  that  a  top-loaded  tube  performs  slightly 
better  than  an  unloaded  tube  (0.59  dB  improvement),  the  top  load  was  omitted  for  mechanical 
reasons. 

Fig.  25  gives  a  schematic  view  of  the  Loran-C  antenna  and  amplifier  assembly.  For 
mechanical  reasons,  the  antenna  and  the  amplifier  were  enclosed  in  a  1/16  in.  thick  fiberglas  tube 
of  5  in.  outer  diameter.  The  antenna  was  constructed  from  a  1/16  in.  thick  aluminum  tube  that  fits 
tightly  against  the  fiberglas  housing.  The  amplifier  was  enclosed  in  a  grounded  aluminum  can  of 
0.965  in.  height  and  the  same  diameter  as  the  antenna.  The  can  included  a  3/16  in.  diameter  hole 
through  its  middle  for  passage  of  the  GPS  antenna  feed  cable.  The  Loran-C  antenna  and  amplifier 
assembly  in  Fig.  25  was  enclosed  in  a  cylindrical  radome  made  of  lexan  (er  =  3)  and  further 

surrounded  by  a  floatation  collar  made  of  syntactic  foam  (er  =  1.5). 
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Since  the  amplifier  required  a  larger  enclosure  than  previously  anticipated  (0.965  in.  height 
instead  of  0.75  in.),  the  antenna  height  was  reduced  to  3.785  in.  while  the  fiberglas  housing 
reduced  the  antenna  diameter  to  approximately  4  7/8  in.  Additionally,  the  tube  was  separated  by 
0.25  in.  from  the  ground  plane.  The  recalculation  of  the  parameters  for  this  smaller  antenna  (h  = 
3.785  in.,  2a  =  4.875  in.,  and  s  =  0.25  in.)  yields  Ca  =  17.8  pF,  heg-=  1.53  in.,  and  Eamp  =  4.60  pV/m. 

From  the  CCIR  tables  [4],  the  atmospheric  noise  (for  a  100  Hz  bandwith)  near  the  coast  of 
Connecticut  ranges  from  a  minimum  of  0.94  nV/m  during  the  Winter  to  a  maximum  of  21  |iV/m 
during  the  Summer.  Thus,  according  to  these  numbers,  the  Loran  antenna/amplifier  system  is 
usually  atmospheric  noise  limited. 

The  effective  height  of  the  Loran  antenna  was  measured  by  comparing  the  received  signal 
strength  from  a  Loran  station  with  that  from  a  VLF  loop  antenna  of  known  effective  height.  Based 
on  this  method,  the  measured  effective  height  was  approximately  1  in.,  smaller  than  the  number 
predicted.  The  smaller  measured  effective  height  is  probably  due  to  the  added  capacitive  loading 
from  the  dielectric  surrounding  the  antenna  and  the  finite  thickness  of  the  connecting  wire. 
Reception  problems  occurred  when  the  floatation  collar  became  wet.  The  measurements  were 
carried  out  over  a  salt  water  pond  in  Fisher’s  Island,  NY. 

SUMMARY  AND  CONCLUSIONS 

This  report  has  provided  a  systematic  theoretical  determination  of  an  optimum  monopole 
antenna  configuration  for  Loran-C  reception  that  fits  within  certain  space  requirements.  The 
antenna  analysis  was  based  on  a  quasi-static  numerical  study  of  a  conducting  body  of  revolution 
above  a  perfectly  conducting  ground  plane.  An  optimum  antenna  is  one  that  results  in  a 
minimum  amplifier  noise  field  as  computed  from  (7). 

An  integral  part  of  this  study  included  a  determination  of  the  amplifier  noise.  On  account  of 
the  small  dimensions  allotted  for  the  Loran  antenna,  the  antenna  capacitance  was  sufficiently 
small  so  that  the  amplifier  was  current  noise  limited  (see  Fig.  4),  resulting  in  a  noise  field  that  is 
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inversely  proportional  to  the  product  of  the  antenna  capacitance  and  effective  height. 
Consequently,  the  study  showed  a  tradeoff  in  the  effort  to  obtain  both  a  large  antenna  capacitance 
and  effective  height.  For  example,  a  structure  which  is  more  tapered  at  the  base  yields  an  increase 
in  effective  height  at  the  expense  of  a  loss  in  capacitance.  The  analysis  showed  that  the  optimum 
monopole  antenna  is  a  hollow  cylinder  of  maximum  allowable  diameter  and  a  small  separation 
from  the  ground  plane.  The  analysis  also  indicated  that  the  addition  of  a  top  load  provides  a  small 
increase  in  both  the  antenna  capacitance  and  effective  height  resulting  in  an  improved 
performance.  However,  because  of  mechanical  constraints,  it  was  decided  to  omit  the  top  load. 

Measurements  performed  on  the  Loran-C  antenna  indicate  a  lower  effective  height  than  the 
predicted  value  because  of  the  finite  thickness  of  the  connecting  wire  and  the  added  shunt 
capacitive  loading  from  the  dielectric  surrounding  the  antenna.  The  computer  program  for  the 
conducting  body  of  revolution  assumes  an  upper-half  space  consisting  of  air  and  does  not  account 
for  any  dielectric  loading.  The  Loran  antenna  performance  degraded  significantly  when  it 
became  wet. 
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Figure  1  Loran-C/GPS  integrated  antenna  system. 
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(b) 


Figure  2.  (  a )  Monopole  antenna  above  a  ground  plane. 

( b )  Thevenin  equivalent  circuit  in  the  receiving  mode. 
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Antenna  Connecting  Amplifier 

Wire  Input 


(b) 


Figure  3.  (  a )  Monopole  antenna  connected  to  an  FET  amplifier, 
(b)  Equivalent  circuit 
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(b) 


Figure  4.  Amplifier  noise  model. 

( a )  Input  noise  voltage  of  the  amplifier. 

( b )  Input  noise  current  In  of  the  amplifier. 
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Figure  5.  Equivalent  circuit  for  determination  of  the  voltage  transfer  function. 
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Figure  6.  Body 


and  coordinate  system. 
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Figure  9.  Conducting  tube  above  a  ground  plane. 
(  a  )  without  top  load 
(  b  )  with  top  load 


Figure  11.  Capacitances  of  several  conducting  tubes  (  h 
ground  plane. 


•  2a  =  3.5  in.  (measured) 
2a  =j4.5_in.  (computed) 
□  2a  =  4.5  in.  (measured) 
2a  =  5.6  in.  (computed) 
K  2a  =  5.6  in.  (measured) 
2a  —  6.2  in.  (computed) 
A  2a  =  8.2  in.  (measured) 


A 


6  in.)  of  different  radii  as  a  function  of  separation  above  a 


Figure  14. 


Capacitances  of  several  top-loaded  conducting  tu 
of  separation  above  a  ground  plane. 


Legend 

2a  =  3.5  in.  (computed) 
•  2a  =  3.5  in.  (measured) 
2a  =  5.6  in.  (computed)^ 
□  2a  =  5.6  in.  (measured) 
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Figure  17.  Top-loaded  conducting  tube  above  a  ground  plane. 
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Figure  20.  Effective  height  of  a  truncated  cone  ( h  =  6  in.,  df  =  5  in.)  above  a  ground  plane  (  s  =  0.1  in.)  as  a  function  of 
the  lower  diameter. 
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Figure  23.  Top-loaded  truncated  cone  above  a  ground  plane. 
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APPENDIX  A 

ANALYSIS  OF  THE  SINGULARITY  IN  THE  INTEGRAND  OF  EQN.  ( 18 ) 


In  this  analysis,  the  following  parametrization  is  valid: 

Z  -  z'  =  (t  -  t’)  COS  Vj 

p  -  p'  =  (t  - 1')  sin  Vj 
where  Vj  is  defined  by  (23).  From  (14), 


•n/2 


da 


sVpp7  I  2  .  2 

/  V  v  +  sin  a 


( A  -  la ) 
(A- lb) 


(A-2) 


where 


y  = 


(p-p')2+(z-z’)2 


1/2 


It-fl 


( A-  3 ) 


2/pp7  2VpfF 

It  is  desired  to  determine  the  above  integral  when  t’  t  (y  — >  0).  As  t'  — >  t,  the  latter  integrand  in 

(A-2)  is  singular  when  a  =  0.  Equation  (A-2)  may  be  expressed  as 

1 


K(u2) 


Ri 


( A-4) 


where 


IrfY)  = 


1-cosa 


A f~2  .  2 

V  Y  +  sin  a 


da 


(A- 5a) 


W  = 


cos  a 


.  2 

V  Y  +  sin  a 


da 


(A- 5b) 


The  integrand  of  Ib(Y)  is  bounded  for  y  =  0  and  a  =  0  while  the  integrand  of  Is(y)  is  singular. 
An  evaluation  of  Ib  at  y  =  0  is  given  as  follows: 
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lim  I^Y)  = 
y-»  0 


rnll 

1  1_cosa  da  -  In 

1 1-cosa' 

1  sin  a 

2 

' 0 

\  sin  a 

a  =  jt/2 

a  =  0 


•  2  \ 

=  lim  In  3in  a-  =  ln(2) 
a  -» 0  U-cosa/ 


The  substitution  of  x  =  sin  a  in  (A-5b)  gives 


(A-  6 ) 


Therefore, 


lim  Uy)  =  -lnjlllj 


Thus  I  has  a  logarithmic  singularity  at  y  =  0.  From  (A-6)  and  (A-7), 


(A-  7 ) 


The  use  of  (A-8)  in  (A-4)  yields 

2 

«-.»  r{  ^y-o 


^-[ln(4)  +  ln(2p)-lnlt-t’l 


(A-8) 


(A— 9) 


A  similar  result  was  obtained  by  Glisson  and  Wilton  [11]  in  their  study  of  electromagnetic 
scattering  from  bodies  of  revolution. 
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APPENDIX  B 
COMPUTER  PROGRAM 

The  double  precision  FORTRAN  computer  program  listed  in  this  appendix  computes  the 
charge  distribution,  input  capacitance,  and  effective  height  of  a  conducting  body  of  revolution 
above  a  conducting  ground  plane.  The  charge  distribution  along  the  body  is  obtained  from  the 
solution  of  (18)  by  the  method  of  moments  as  described  in  this  report.  The  capacitance  is  given  by 
(19)  while  the  effective  height  is  obtained  from  (46).  In  addition,  the  amplifier  noise  field  is 
determined  from  (7). 

The  endpoints  of  the  segments  approximating  the  generating  curve  ( (pj,^),  i  =  0,  1,  2, ...,  N)  of 
the  body  of  revolution  are  read  from  a  data  file.  The  coefficient  matrix  S  has  elements  given  by 
(28)  for  i  *  j  and  (34)  for  i  =  j  that  involve  integrals  with  the  following  integrands: 

FINT  =  K(U|j)  ,  i*j  (B- 1) 


FINTD  = 


( B-2) 


K(wf:) 

FINTI  =  jl.  ( B-3) 

RU 

where  all  variables  were  previously  defined  while  K  is  the  complete  elliptic  integral  of  the  first 
kind.  The  function  ELIPT  computes  K  by  a  polynomial  approximation  given  in  [5].  The  above 
integrands  are  functions  of  x'.  Special  care  is  taken  so  that  FINTD  is  not  evaluated  at  x'  =  Aj/2. 

The  integrals  are  computed  by  calling  the  function  SIMP.  SIMP  computes  the  integral  of  a  real 
function  with  finite  limits  by  Simpson's  rule  with  automatic  halving.  The  relative  error  for  the 
integrations  is  given  by  REL.  In  this  study,  REL  was  set  at  1CT4  for  4-place  accuracy. 

The  linear  system  of  equations,  given  in  matrix  form  by  (35),  is  solved  for  the  unknown 
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charge  density  vector  Q  by  calling  the  subroutines  FACTR  and  SOLVE.  FACTR  factors  the 
coefficient  matrix  S  into  the  product  of  a  lower  triangular  matrix  and  an  upper  triangular  matrix 
using  the  Gauss-Doolittle  algorithm  [12].  The  factored  matrix  is  used  by  SOLVE  to  determine  the 
solution  of  the  matrix  equation.  The  excitation  vector  is  read  into  Q  and  is  overwritten  with  the 
charge  density  vector  after  calling  SOLVE.  FACTR  and  SOLVE  are  part  of  the  Numerical 
Electromagnetics  Code  [13]. 
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0  ***************************************************************  *  * *********** 
C  THIS  PROGRAM  COMPUTES  THE  ELECTROSTATIC  CHARGE  PER  UNIT  LENGTH  OF  A 
C  CONDUCTING  BODY  OF  REVOLUTION  CHARGED  TO  A  CONSTANT  POTENTIAL  AND  LYING 
C  ABOVE  A  CONDUCTING  GROUND  PLANE.  AN  INTEGRAL  EQUATION  IS  SOLVED  FOR  THE 
C  CHARGE  DENSITY  BY  THE  METHOD  OF  MOMENTS  USING  PULSE  EXPANSION  FUNCTIONS 
C  AND  POINT  MATCHING.  THE  RESULT  IS  APPLIED  TO  EVALUATE  THE  CAPACITANCE, 

C  EFFECTIVE  HEIGHT,  AND  AMPLIFIER  NOISE  FIELD  OF  AN  ELECTRICALLY  SMALL 
C  ANTENNA  CONSISTING  OF  THE  BODY  OF  REVOLUTION  AND  AN  INFINITESIMALLY  THIN 
C  WIRE  CONNECTING  THE  BODY  TO  A  RECEIVER  LOCATED  BELOW  THE  GROUND  PLANE. 

C 

C  THE  INPUT  VARIABLES  ARE  GIVEN  AS  FOLLOWS: 

C  REL  =  RELATIVE  ERROR  FOR  INTEGRATIONS 

C  VOLT  =  POTENTIAL  OF  BODY  IN  VOLTS 

C  CW  =■  INTERCONNECTION  CAPACITANCE  IN  PICOFARADS 

C  NVOLT  =  AMPLIFIER  NOISE  VOLTAGE  IN  NANOVOLTS 

C  NCUR  =  AMPLIFIER  NOISE  CURRENT  IN  PICOAMPS 

C  N  =  NO.  OF  BASIS  FUNCTIONS 

C  RHOB(I)  =  RADIAL  COORDINATE  OF  ENDPOINT  I  IN  INCHES 
C  ZB (I)  =  AXIAL  COORDINATE  OF  ENDPOINT  I  IN  INCHES 

C 

C  INTERMEDIATE  VARIABLES  ARE  GIVEN  AS  FOLLOWS: 

C  RHO(I)  =  RADIAL  COORDINATE  OF  MIDPOINT  OF  SEGMENT  I  IN  METERS 
C  Z(I)  =  AXIAL  COORDINATE  OF  MIDPOINT  OF  SEGMENT  I  IN  METERS 

C  T (I)  =  T  COORDINATE  OF  MIDPOINT  OF  SEGMENT  I  IN  METERS 

C  D(I)  »  LENGTH  OF  SEGMENT  I  IN  METERS 

C  V (I)  =  ANGLE  BETWEEN  SEGMENT  I  AND  THE  Z-AXIS 

C  S(I,J)  =  <I,J)  ELEMENT  OF  COEFFICIENT  MATRIX 

C 

C  THE  OUTPUT  VARIABLES  ARE  GIVEN  AS  FOLLOWS: 

C  Q (I)  =  CHARGE  DENSITY  ALONG  MIDPOINT  OF  SEGMENT  I  IN  COUL/M 

C  CAP  -  ANTENNA  CAPACITANCE  IN  FARADS 

C  EH  =  ANTENNA  EFFECTIVE  HEIGHT  IN  METERS 

C  EAMP  =  AMPLIFIER  NOISE  FIELD  IN  VOLTS/M 

Q  **************************************************************************** 

c 

IMPLICIT  REAL* 8  (A-H,0-Z) 

REAL* 8  NCUR, NVOLT 

DIMENSION  RHO (200) , RHOB (0:200) , Z (200) ,  ZB (0 : 200 )  ,T(200)  ,V(200)  ,D (200) 
DIMENSION  Q (200) ,S (200,200) , IP (200) 

EXTERNAL  FINT, FINTD, FINTI , ELIPT 
COMMON  /Cl/  RI , RJ1 , ZI, ZJ1 , CVJ, SVJ,DI2 
C 

C  ****  CONSTANTS  **** 

C 

PI=3. 1415926535897 9D0 

E0=8 . 854D-12  !  PERMITTIVITY  OF  FREE  SPACE 

VOLT=l .DO  !  POTENTIAL  IN  VOLTS 
FR=l.D+5  !  FREQUENCY 
C 

C  ****  READ  INPUT  DATA  **** 

C 

WRI TE ( 6 , 1 ) 

WRITE (7,1) 

1  FORMAT (1 OX, 'BODY  OF  REVOLUTION  ABOVE  A  CONDUCTING  GROUND  PLANE') 

WRI TE ( 6 , 1 0 ) 

10  FORMAT (/,5X, 'ENTER  RELATIVE  ERROR  FOR  INTEGRATIONS  :  ' , $) 

READ ( 6 , * )  REL 
WRITE (7, 15)  REL 
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15 

20 


25 

30 


35 


40 


45 


50 


55 


60 

65 


80 

70 


100 


110 


FORMAT (1 OX, 'RELATIVE.  ERROR  FOR  INTEGRATIONS  =  ?,E12.5) 

WRITE (6, 20) 

FORMAT  (/,  5X, 'ENTER  POTENTIAL  OF  BODY  IN  VOLTS  :  ' , $) 

READ  ( 6  ,  * )  VOLT 
WRITE (7, 25)  VOLT 

FORMAT </,5X, 'POTENTIAL  OF  BODY  =  ',E12.5, '  VOLTS') 

WRITE (6, 30) 

FORMAT </,5X, 'ENTER  INTERCONNECTION  CAPACITANCE  IN  PICOFARADS', 

*  '  :  ',$) 

READ  ( 6  ,  * )  CW 
WRITE (7, 35)  CW 

FORMAT (/, 5X, ' INTERCONNECTION  CAPACITANCE  =  ' , E12 . 5, '  PICOFARADS ' ) 
CW=CW*1.D-12 


WRITE (6, 40) 

FORMAT (/, 5X, 'ENTER  AMPLIFIER 

READ (6, *)  NVOLT 

WRITE (7, 45)  NVOLT 

FORMAT (/, 5X, 'AMPLIFIER  NOISE 

NVOLT=NVOLT*l .D— 9 

WRITE (6, 50) 

FORMAT (/, 5X, 'ENTER  AMPLIFIER 

READ (6,*)  NCUR 

WRITE (7, 55)  NCUR 

FORMAT (/,5X, 'AMPLIFIER  NOISE 

NCUR=NCUR* 1 . D-12 

READ  (8,60)  N  !  READ  NO.  OF 

FORMAT (IX, 13) 

WRITE (7, 65) 

FORMAT (/, 5X, 'NO.  OF  SEGMENTS 
DO  70  1=0, N  !  READ  SEGMENT 
READ (8,80)  RHOB(I) ,ZB(I) 
FORMAT (IX, E22 . 15, 5X, E22 . 15) 
CONTINUE 


NOISE  VOLTAGE  IN  NANOVOLTS  :  ' , $) 

VOLTAGE  =  ' , E12 . 5, '  NANOVOLTS ' ) 

NOISE  CURRENT  IN  PICOAMPS  :  ',$) 

CURRENT  =  ' , E12 . 5, '  PICOAMPS') 
SEGMENTS 

-  ’,13) 

ENDPOINTS 


DO  100  1=0, N  !  CONVERT  SEGMENT  ENDPOINTS  TO  METERS 
RHOB (I ) =RHOB (I) *2 . 54D-2 
ZB (I ) =ZB  (I ) *2 . 54D-2 
DO  110  1=1, N 

RHO (I) = (RHOB (1-1) +RHOB (I) )*.5D0 
Z (I) = ( ZB (1-1) +ZB (I) ) *.5D0 

V (I) =ATAN2 (RHOB (I) -RHOB (1-1) , ZB (I) -ZB (1-1) ) 
D(I)=SQRT( (RHOB (I) -RHOB (1-1) ) **2+ (ZB (I) -ZB (I— 1) ) **2) 
IF  (I.EQ.l)  THEN 
T (I ) = . 5D0*D (I) 

ELSE 

T (I ) =T (1-1) + . 5D0* (D (1-1) +D (I) ) 

END  IF 
CONTINUE 


****  EVALUATE  COEFFICIENT  MATRIX  **** 


DO  200  1=1, N 
RI=RHO ( I ) 

ZI=Z (I) 

DI2=. 5D0*D  (I) 
DO  210  J=1,N 
RJl=RHOB (J-l) 
ZJ1=ZB( J-l) 
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CVJ=COS (V(J)  ) 

SVJ=SIN (V ( J)  ) 

SIMAGE=SIMP (FINTI ,  0  .  DO ,  D  ( J) ,REL) 

IF  (I.EQ.J)  GO  TO  220 

S (I,  J) =SIMP (FINT, 0 .DO, D ( J) , REL) -SIMAGE 
GO  TO  210 
220  EP=Dl2*l.D-2 
EPL=DI2-EP 
EPU=DI2+EP 

ST-SIMP (FINTD, 0 .DO, EPL, REL) +SIMP (FINTD, EPU, D ( I )  ,  REL) 

DO  230  L=l, 5 
EP=EP* . 1D0 
EPL1=DI2-EP 
EPU1=DI2+EP 

ST1=SIMP (FINTD, EPL, EPL1, REL) +SIMP (FINTD, EPU1 , EPU, REL) +ST 
IF  (ABS (ST1-ST) .LT.REL*ABS (ST) )  THEN 

S (1,1) =ST1+DI2* (1. DO-LOG (DI2) ) /RHO (I ) -SIMAGE 
GO  TO  210 
ELSE 

ST=ST1 
EPL=EPL1 
EPU=EPU1 
END  IF 

230  CONTINUE 

S (1,1) =ST1+DI2* (1 .DO -LOG (DI2) ) /RHO(I) -SIMAGE 
WRITE (6,240)  I, I 

240  FORMAT (//, 10X, ' S ( ' , 13, ', ',13, ' )  I S  NOT  ACCURATE !' ) 

210  CONTINUE 
200  CONTINUE 

****  COMPUTE  INPUT  EXCITATION  VECTOR  **** 

TEMP=2.D0*PI*PI*E0*VOLT 
DO  250  1=1, N 
Q (I) =TEMP 

****  COMPUTE  CHARGE  DISTRIBUTION  **** 

CALL  FACTR (N, S, IP, 200) 

CALL  SOLVE (N,S, IP, Q, 200) 

****  COMPUTE  CAPACITANCE  **** 

QT=0 . DO 
DO  260  1=1, N 

QT=QT+Q (I) *D (I)  !  TOTAL  CHARGE  ALONG  BODY 

CAP=QT /VOLT 

****  COMPUTE  EFFECTIVE  HEIGHT  **** 

EH=0.D0 
DO  270  1=1, N 
270  EH=EH+Q { I ) *D ( I ) * ( Z ( I )  -ZB  ( 0 )  ) 

EH=EH/QT 

EH=EH+ZB ( 0 )  !  ACCOUNT  FOR  HEIGHT  OF  CONNECTING  WIRE 

C 

EAMP= (NVOLT* (1 .D0+CW/CAP) ) **2+(NCUR/ (2 .D0*PI*FR*CAP) ) **2 
EAMP=SQRT (EAMP) /EH  !  AMPLIFIER  NOISE  FIELD  IN  VOLTS/METER 
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****  PRINT  OUTPUT  **** 


300 


310 

320 

330 


WRITE (7,300)  (I, T (I) /2 . 54D-2, 1,Q (I) ,1=1, N) 

FORMAT (IX, 'T(', 13, ')  =  ',E14.7,»  in. ' , 5X, ' Q ( • f 13, ' )  =  ',E14.7 
'COUL/M* ) 

WRITE (7, 310)  CAP*1.D+12 

FORMAT <10X, 'CAPACITANCE  =  ' ,E14.7,'  pF') 

WRITE (7, 320)  EH/2.54D-2 

FORMAT (10X, 'EFFECTIVE  HEIGHT  =  ' ,E14.7,'  in.') 

WRITE (7, 330)  EAMP*1 . D+6 

FORMAT (10X, 'AMPLIFIER  NOISE  FIELD  =  *,E14.7,'  microvolts /m' ) 

STOP 

END 


FUNCTION  FINT(T) 

IMPLICIT  REAL* 8  (A-H,M-Z) 

COMMON  /Cl/  RI,RJ1,ZI,ZJ1,CVJ,SVJ,DI2 
R2=SQRT { (RI+RJ1+T*SVJ) **2+ (ZI-ZJ1-T*CVJ) **2) 
U2=4.D0*RI* (RJ1+T*SVJ) / (R2**2) 

M1=1.D0-U2 
FINT=ELIPT (Ml) /R2 
RETURN 
END 


FUNCTION  FINTD(T) 

IMPLICIT  REAL* 8  (A-H,M-Z) 

COMMON  /Cl/  RI , RJ1 , ZI ,  Z J1 , CVJ, SVJ, DI2 
R2=SQRT ( <RI+RJ1+T*SVJ) **2+ (ZI-ZJ1-T*CVJ) **2) 
U2=4.D0*RI* (RJ1+T*SVJ)  / (R2**2) 

Ml=l . DO— U2 

FINTD=ELIPT (Ml) /R2+. 5DO*LOG (ABS (DI2-T) ) /RI 

RETURN 

END 


FUNCTION  FINTI(T) 

IMPLICIT  REAL* 8  (A-H,M-Z) 

COMMON  /Cl/  RI,RJ1,ZI,ZJ1,CVJ, SVJ,DI2 
R3=SQRT ( (RI+RJ1+T*SVJ) **2+ (ZI+Z J1+T*CVJ) **2) 
G2=4 .D0*RI* (RJ1+T*SVJ) / (R3**2) 

Ml=l .D0-G2 
FINTI=ELIPT (Ml) /R3 
RETURN 
END 


FUNCTION  ELIPT (AMI ) 

c  ******-****************************************************** 

C  COMPLETE  ELLIPTIC  INTEGRAL  OF  THE  FIRST  KIND, K (M) , AS  DEFINED 

C  IN  THE  REFERENCE  BELOW, WHERE  AM1=1-M  . 

C 

C  REFERENCE:  HANDBOOK  OF  MATHEMATICAL  FUNCTIONS  (AMS-55) 

C  AUTHORS:  ABRAMOWITZ  AND  STEGUN 

C  EQ.  17.3.34 
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IMPLICIT  REAL* 8  (A-H,0-Z) 

DATA  AO,  Al,  A2,  A3,  A4,  BO,  Bl,  B2,  B3 , B4/1 . 38629436112D0, 

*  .9666344259D-1, . 3590092383D-1, . 3742563713D-1, 

*  . 1451196212D-1, . 5D0, . 12498593597D0, . 6880248576D-1 , 

*  . 3328355346D-1, . 441787012D-2/ 

A=A0+A1  *  AMI 

B=B0+B1*AM1 

IF  (AM1.LT.1.D-18)  GO  TO  10 

AM12=AM1*AM1 

A=A+A2  * AMI 2 

B=B+B2  * AMI 2 

IF  (AM1.LT.1.D-12)  GO  TO  10 

AM13=AM12*AM1 

A=A+ A3  * AMI  3 

B=B+B3*AM13 

IF  (AM1.LT.1.D-9)  GO  TO  10 
AM14=AM13*AM1 
A=A+A4*AM14 
B=B+B4*AM14 
10  CONTINUE 

ELIPT=A-B*  LOG (AMI ) 

RETURN 

END 


FUNCTION  SIMP (F, A,B,REL) 

**************************************************** 

INTEGRATION  BY  SIMPSON'S  RULE  WITH  AUTOMATIC  HALVING 
**************************************************** 

F  =  REAL  FUNCTION  TO  BE  INTEGRATED 
A  =  LOWER  LIMIT  OF  INTEGRATION 
B  =  UPPER  LIMIT  OF  INTEGRATION 
REL  =  RELATIVE  ERROR  INDICATOR 

IMPLICIT  REAL* 8  <A-H,0-Z) 

ST=0 -DO 
SA=F (A) 

SB=F  (B) 

N=2 

H= (B-A) * . 5D0 
S=  (SA+SB)  * .  5D0 
DO  10  1=1,7000 
T=0.D0 

DO  20  J=1 , N-l , 2 
20  T=T+F (A+H* J) 

S=S+T  1 

V=(S4T) *H* . 666666666 6 66 6667D0 

IF  (ABS (V-ST) .LT.REL*ABS (ST) )  GO  TO  30 

ST=V 

N=N4N 

10  H=H* . 5D0 
30  SIMP=V 
RETURN 
END 
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SUBROUTINE  FACTR (N, A, IP,NDIM) 

Q  ****************************************************************** 

C  SUBROUTINE  TO  FACTOR  A  MATRIX  INTO  A  UNIT  LOWER  TRIANGULAR  MATRIX 

C  AND  AN  UPPER  TRIANGULAR  MATRIX  USING  THE  GAUSS-DOOLITTLE  ALGORITHM 

C  PRESENTED  ON  PAGES  411-416  OF  A.  RALSTON — A  FIRST  COURSE  IN 

C  NUMERICAL  ANALYSIS.  COMMENTS  BELOW  REFER  TO  COMMENTS  IN  RALSTONS 

C  TEXT.  (MATRIX  TRANSPOSED. 

Q  ****************************************************************** 

IMPLICIT  REAL* 8  (A-H,0-Z) 

DIMENSION  A (NDIM,NDIM) ,  IP (NDIM) 

COMMON  /SCRATM/  D(600) 

INTEGER  R, RM1 , RP1 , P J, PR 
IFLG=0 
DO  9  R=1 , N 

STEP  1 

DO  1  K=1,N 
D(K)=A(R,K) 

CONTINUE 

STEPS  2  AND  3 

RM1=R— 1 

IF  (RM1.LT.1)  GO  TO  4 
DO  3  J=1 , RM1 
PJ=IP(J) 

ARJ=D (P J) 

A (R, J) =ARJ 
D(PJ)=D(J) 

JP1-J+1 
DO  2  I=JP1,N 
D (I) =D (I) -A (J, I) *ARJ 
CONTINUE 
CONTINUE 
CONTINUE 

STEP  4 


DMAX=D (R) **2 
IP (R) =R 
RP1-R+1 

IF  (RP1.GT.N)  GO  TO  6 
DO  5  I=RP1, N 
ELMAG=D ( I ) *  *  2 

IF  (ELMAG.LT.DMAX)  GO  TO  5 

DMAX=ELMAG 

IP (R)=I 

CONTINUE 

CONTINUE 

IF  (DMAX.LT.l.D-10)  IFLG=1 
PR=IP (R) 

A  (R,  R)  =D  (PR) 

D (PR) =D (R) 

STEP  5 

IF  (RP1.GT.N)  GO  TO  8 
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ARJ=1 .  DO/A  (R,  R) 

DO  7  I=RP1,N 
A(R,  I)  =D  (I)  *ARJ 

7  CONTINUE 

8  CONTINUE 

IF  (IFLG.EQ.O)  GO  TO  9 
PRINT  10,  R, DMAX 
IFLG=0 

9  CONTINUE 
RETURN 

FORMAT  (1H  , 6HPIVOT (,I3,2H)=, El 6 . 8 ) 
END 


SUBROUTINE  SOLVE  (N, A, IP, B,NDIM) 

**************************************************************** 
SUBROUTINE  TO  SOLVE  THE  MATRIX  EQUATION  LU*X=B  WHERE  L  IS  A  UNIT 
LOWER  TRIANGULAR  MATRIX  AND  U  IS  AN  UPPER  TRIANGULAR  MATRIX  BOTH 
OF  WHICH  ARE  STORED  IN  A.  THE  RHS  VECTOR  B  IS  INPUT  AND  THE 
SOLUTION  IS  RETURNED  THROUGH  VECTOR  B.  (MATRIX  TRANSPOSED. 

**************************************************************** 

IMPLICIT  REAL* 8  (A-H,0-Z) 

INTEGER  PI 

COMMON  /SCRATM/  Y(600) 

DIMENSION  A (NDIM,NDIM) ,  IP (NDIM) ,  B (NDIM) 

FORWARD  SUBSTITUTION 

DO  3  1=1, N 
PI=IP(I) 

Y(I)=B(PI) 

B(PI)=B(I) 

IP1=I+1 

IF  (IP1.GT.N)  GO  TO  2 
DO  1  J=IP1,N 
B  ( J)  =B  ( J)  -A  (I,  J)  *Y  (I ) 

CONTINUE 
CONTINUE 
CONTINUE 

BACKWARD  SUBSTITUTION 

DO  6  K=1,N 
I=N— K+l 
SUM=0 . DO 
IP1=I+1 

IF  (IP1.GT.N)  GO  TO  5 
DO  4  J=IP1,N 
SUM=SUM+A ( J, I ) *B ( J) 

4  CONTINUE 

5  CONTINUE 

B(I)  =  (Y(I)  -SUM)  /A  ( I ,  I ) 

6  CONTINUE 
RETURN 
END 
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APPENDIX  C 

APPROXIMATE  EXPRESSION  FOR  THE  CAPACITANCE  OF  A 
SMALL  TUBULAR  ANTENNA 

Consider  the  tubular  dipole  shown  in  Fig.  C-l(a).  It  consists  of  two  perfectly  conducting  hollow 
tubes  each  of  length  h  and  radius  a  and  separated  by  a  distance  2s.  The  charge  distribution  along 
the  dipole  is  rotationally  symmetric.  To  compute  the  capacitance  of  the  tubular  dipole,  consider  the 
coplanar  stripline  geometry  in  Fig.  C-l(b).  It  consists  of  two  strips  each  of  width  h  and  length  2na 
and  separated  by  a  distance  2s.  The  capacitance  per  unit  length  CL  of  the  coplanar  stripline  is 

given  by  [14] 


c  _eJ^2 

~  eo  K(k) 


where 


k  = 


s  +  h 


(C-l) 


(C-2) 


K(k)  is  the  complete  elliptic  integral  of  the  first  kind  [5]  and  k’  =  (1  -  k2)172.  The  capacitance  of  the 
stripline  in  Fig.  C-l(b)  is  obtained  by  multiplying  CL  by  the  length  2na.  Therefore,  by 

approximating  the  tubular  dipole  in  Fig.  C-l(a)  by  the  coplanar  stripline  in  Fig.  C-l(b),  the  dipole 
capacitance  is 


c  =  27taeoWT  (C_3) 

The  capacitance  of  a  tubular  monopole  of  height  h  and  radius  a  and  separated  by  a  distance  s  above 
an  infinite  ground  plane  is  given  by  twice  the  value  in  (C-3). 
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Figure  C  -  1  (a)  Tubular  dipole. 

(b)  Equivalent  coplanar  stripline. 
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